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and Power Allocation Algorithms in a Relay with 
Partial Channel State Information 
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Abstract —The performance of a dual-hop MIMO relay net¬ 
work is studied in this paper. The relay is assumed to have access 
to the statistical channel state information of its preceding and 
following channels and it is assumed that fading at the antennas 
of the relay is correlated. The cumulative density function (cdf) 
of the received SNR at the destination is first studied and 
closed-form expressions are derived for the asymptotic cases of 
the fully-correlated and non-correlated scenarios; moreover, the 
statistical characteristics of the SNR are further studied and an 
approximate cdf of the SNR is derived for arbitrary correlation. 
The cdf is a multipartite function which does not easily lend 
itself to further mathematical calculations, e.g., rate optimization. 
However, we use it to propose a simple power allocation algorithm 
which we call “proportional power allocation”. The algorithm 
is explained in detail for the case of two antennas and three 
antennas at the relay and the extension of the algorithm to a 
relay with an arbitrary number of the antennas is discussed. 
Although the proposed method is not claimed to be optimal, the 
result is indistinguishable from the benchmark obtained using 
exhaustive search. The simplicity of the algorithm combined with 
its precision is indeed attractive from the practical point of view. 


Index Terms —Cooperative Communications, AmpUfy-and- 
Forward, Multi-Antenna Relay, Statistical Channel State Infor¬ 
mation, Largest Eigenmode Relaying 

I. Introduction 

I T was proved by van der Meulen in |[T| that the capacity 
of a three node communication system can, potentially, 
be larger than the capacity of a point-to-point communication 
system. Consequently, analysis of communication systems in 
which the transceiver nodes cooperatively transmit their data 
to an intended final receiver has been a rather active field 
of research in last decade and numerous papers including, 
e.g., 0-0 have investigated cooperative communication 
systems. During the infancy of the concept of cooperative 
communications, substantial work was carried out, investigat¬ 
ing cooperating nodes with single antennas. Several promising 
relaying protocols were proposed; among them. Amplify and 
Forward (AF) is intensively studied in the literature; hence, in 
this paper, we will focus on AF relaying too. 

Employing multiple antennas in communication nodes is 
another technique proved to be capable of enhancing trans¬ 
mission rates (see e.g., 0). Employing multiple antennas in 
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the nodes of a cooperative communication system has been 
an active research trend during the last few years. Assuming 
multiple antennas at the relay, one major task is to design a 
suitable amplihcation matrix in the relay. Indeed, depending on 
the available Channel State Information (CSI) at the relay, the 
amplification matrix can, potentially, be different. Moreover 
different communication systems can demand the optimization 
of different desired performance measure; hence, different 
“optimal” relaying protocols will exist; for instance, the non- 
regenerative relaying matrices, e.g., in fTl-pOl, are designed 
to minimize Mean Square Error (MSE) but other relaying 
matrices, e.g., in im -p3|, are assumed to maximize the 
achievable rates. 


Assuming statistical CSI at a transceiver node is interest¬ 
ing from a practical point of view; in particular, in rapidly 
changing channels, assuming perfect CSI in a relay node 
is, indeed, unrealistic, hence, a large body of the literature 
investigates AE cooperative systems wherein a single antenna 
relay node has access only to the statistical CSI (see, e.g., 
fT^). Note that single antenna AE relaying systems, with 
statistical CSI at the relay, are usually referred to as “fixed 
gain” AE relaying. In spite of the importance of cooperative 
communication systems with statistical CSI knowledge, very 
few papers consider the problem when the relay node is 
equipped with multiple antennas. Moreover, except pT| , we 
are not aware of any other paper assuming fading correlation 
in the relay when only the covariance of the channels is known 
to the relay. Note that fading correlation at a transceiver can 
be due to an unobstructed node or space limits at the node 
which forces the antennas to be closely located. Justifications 
to assume transceivers with fading correlation can be found in 

im, (ID. 

The first contribution of this paper is to provide a statistical 
analysis of the received Signal to Noise Ratio (SNR) at the 
destination. There are two major motivations for studying the 
statistical characteristics of the SNR: 


• Outage probability is directly related to received SNR. 
Indeed, the cumulative density function (cdf) of the SNR 
corresponds to the outage probability, and so, the cdf of 
SNR will be derived in this paper. 

• By deriving the cdf of SNR, the mathematical complexity 
of direct maximization of the achievable rate (i.e., “op¬ 
timal” power allocation) will be revealed. It will be an 
excellent motivation for devising alternative approaches 
with reasonable complexity. 
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Accordingly, the second major contribution of the paper is 
to study the problem of power allocation in the relay and 
hence to devise a new and simple power allocation scheme for 
multi-antenna relays. (11), in)’iiD’i[i3 consider the similar 
problem of the power allocation in the relay when statistical 
CSI is available in the nodes (either the source or relay 
nodes). However, whi le fTS) considers the high SNR regime 
of the system, |14|, |19| assumes correlation at the source 
node. In we study a cooperative communication system 
wherein the relay node is equipped with multiple antennas 
that are spatially correlated. The considered system is studied 
only at low SNR and it is proved that Largest Eigenmode 
Relaying (LER) is the optimal transmission method at low 
SNR; however, the system is not studied in the moderate and 
high SNR region. To the best of our knowledge, the design of 
an amplihcation matrix in an AE cooperative system where 
only the statistical CSI is known to the relay is an open 
problem, and one that will be tackled in this paper. We provide 
a scheme which operates in the regime beyond that where LER 
is optimal, and whose performance is indistinguishable from 
the benchmark provided by exhaustive search. 

This paper is organized as follows: In Section]^ the system 
model is introduced and some preliminary existing results are 
recalled. Section III deals with characterizing the statistics of 
the SNR at the destination. In Section IV a simple power 
allocation algorithm is introduced for a relay with only two 
antennas; the proposed algorithm is called “proportional power 
allocation” and has been extended for a system with multiple 
antenna relay node in section IV] and VI and, hnally, the results 
are summarized in Section IVIIl 


IT System Model AND Preliminaries 
A. Notation 

Matrices are represented by boldface upper cases (H). Col¬ 
umn and row vectors are denoted by boldface lower cases (h), 
and hi indicates the i-th element of h. The superscript 
stands for Hermitian transposition. We refer to the identity 
matrix by I. The expectation operation is indicated by E{-}, 
the probability of a random variable is indicated by P(-) and 
fx{x) is reserved for probability density functions (pdf) of 
random variable X\ As represents a diagonal matrix with 
elements organized in descending order and Af denotes the i- 
th diagonal element of A^. Eor simplicity of notation, (Af 
is abbreviated by Af^. The trace of a matrix is denoted by 
Tr(-). 


B. System Model 

In this paper a dual hop, half duplex MIMO communication 
system is investigated. Assume a source node (equipped with 
ns antennas) transmits data to a single antenna destination 
via an intermediate relay node which has nR antennas. The 
proposed system models the downlink of a wide range of 
communication systems in which the user terminal is equipped 
with single antenna due to space limitation, for instance, 
cellular networks or sensor networks. Moreover, fixed-gain 
AE cooperative systems with multiple antennas at the relay 
is an open problem which has received little attention and so 


the proposed system model is a good step forward for under¬ 
standing hxed gain AE systems. It is assumed that a direct 
link between the source and the destination is not available. 
The half duplex constraint is accomplished by time sharing 
between the source and the relay; i.e. each transmission period 
is divided into two time slots: the source transmits during the 
hrst time slot and the relay during the second one. The relay 
remains silent during the source transmission and vice versa. 
It is assumed that the source does not have access to any 
statistical or instantaneous channel state information (CSI). 
Moreover, it is assumed that the antennas in the source node 
are sufficiently far apart and so no correlation is assumed at 
the source . The signal received at the relay iy^ due to the 
source transmission is given by 

y^ = HiX + WYi ( 1 ) 

where the ttr x ns matrix H i represents the channel between 
the source and the relay. With Ps the power constraint of the 
source, the column vector x is the signal transmitted from 
the source with Q — 'E{xx^) = ^Ins the column vector 
Wr represents the receiver noise in the relay with elements 
independently drawn from a complex Gaussian random vari¬ 
able with variance Nq. In this paper, it is assumed that spatial 
correlation occurs at the relay; the correlation can be due to 
space limit in the relay or due to fading correlation due to 
unobstructed relay node. S represents the correlation matrix 
at the relay and therefore, using the Kronecker model, H i can 
be written as 

= (2) 

where elements of Hi^j are i.i.d., zero mean, unit variance 
complex Gaussian random variables, independent of each 
other. The relay multiplies j/r by the gain matrix F and 
forwards it to the destination. Then, the received signal at the 
destination is 

2/d = h2FyR + WB (3) 

= h2FHiXh2FwR +w-D 

where the row vector h ,2 indicates the channel between the re¬ 
lay and the destination; wd represents the receiver noise at the 
destination. Eor simplicity, we assume that wjj and iur are sta¬ 
tistically independent and identical, i.e. = ^o,wr = ^o- 

Due to the spatial correlation S at the relay, one can factorize 
h .2 as 

h.2 = h2w'^^ (4) 

where elements of h 2 w are i.i.d., zero mean, unit variance 
complex Gaussian random variables, independent of each 
other. Justihcation to assume transceivers with spatial correla¬ 
tion can be found in mg, (U). The correlation matrix S in 
the relay is decomposed using spectral decomposition as 

S = (5) 

where is a unitary matrix whose columns are 

the eigenvectors corresponding to S, and As is a 
diagonal matrix with the eigenvalues of S in de¬ 
creasing order, i.e. As = diag[Af, Af, • • • , Af^^] where 
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> Af > • • • > > 0. Moreover, some of Af s can pos¬ 

sibly be zero. 


C. Preliminaries 


Ergodic capacity is one of the main performance criterion 
investigated in this paper. Using ([^, the ergodic capacity of 
the system is defined as 

Cav = l max E{C{Hi,h 2 ,F)} (6) 

F:E{||Fy^||2}<PR 


where C{Hi, h 2 , F) is the conditional transmission rate. For 
simplicity of notation, C{Hi,h 2 , F) is abbreviated by C{-) 
in the rest of the paper. Assuming perfect CSI of i?i and H 2 
at the destination and Q = ^Ins (equal transmit power from 
each antenna in the source, because no channel knowledge is 
available there), the conditional mutual information C'( ) for 
given channel matrices is 


C(-) = log(l + 


Ps h 2 FHiH^F^h^ 
ns No{l + h2FF^h^) 


(7) 


where A) 3 (l -|- h 2 FF^h^) is the total equivalent noise power 
which is assumed to remain constant for coherence time: we 
make a block fading assumption. In pT) , F is found to be 
symmetric as 

F = G5 (8) 


where the gain matrix G is derived as 

G = U^AgU^ (9) 

where Ue is the unitary matrix dehned in 0 and 
Ag = diag[Af, A^, • • • , A^^^]. Note that Af values are to be 
specihed according to the power constraint of the relay so that 
the maximization in (|^ is accomplished; indeed, this is one 
of the main tasks to be handled in this paper. 

Assuming 0 , 0,0 and 0 - the power constraint in the 
relay (i.e. E{|jG5yR|p} < Pr in is 

PsTr(AEAG) + iVoTr(AG) = Pr (10) 


Note that the capacity will be achieved by consuming the 
entire power at the relay, and so, we assume -equality- in ( [T0| ) 
instead of -inequality-. By combining 0, 0, 0, 0 and 
assuming 7 s = Ps/Aq, one can write 0 as follows 


C(-) = log(l + 


7s/12u)AeA, 




ns(l + h 2 wAY,AGh 2 ^) 


) ( 11 ) 


7d 


where 70 represents received SNR in the destination which is 
a function of h 2 w, the correlation eigenvalues matrix 

Ae and the eigenvalues of the G matrix, i.e. Aq. 70 can be 


*The results in this subsection are taken from |ll| . Interested reader is 
recommended to read (0. for a complete proof of the the ideas and the 
derivations; however, in order to make the paper self-contained and also for 
consistency of the notation, the relevant results from m are provided in this 
subsection. 


simplihed according to 


7d = 


7s ^ I h2wAQA'^h\^ i I 

i=l 

ns(l + h 2 wAGA-sh 2 ^) 


(12) 


i-2w^^G-^^Y.i>'2w) 

where hi^^i represents the ith column of ffiu,. Let us assume 


Xj — \h2w,j\ 
1 

ns “ 


(13a) 

(13b) 


It is proved in im Appendix 1] that 7 d in ([0 can be 
further simplihed to 


E 


7d = 7sV X 




1 + E 


(14) 




where n is the minimum of ur and Number of Non-Zero (nnz) 



K = min(nR,nnz(A^)) 


(15) 


Note that Y and Xj correspond to the S-R and R-D channels, 
respectively. The random variable 


E AfAf E 

X = -i-T- „(,) 

1 + E Af Afx, 

i=i 

in ([0 incorporates the effect of the R-D link as well as the 
effect of power allocation due to Xf. Furthermore, since we 
assume Rayleigh fading in both the S-R and R-D links, hence, 
Xj is exponentially distributed with unit mean and Y has an 
Erlang-distribution with rate and shape parameters equal to 
ns- 

Although it is proposed in E) that the optimal G should be 
diagonalized according to 0 where Ag is a diagonal matrix 
with its components organized in descending order, an optimal 
power allocation method to distribute relay’s transmit power 
among different A^s is not discussed. That is still an open 
problem but will be addressed in this paper. 


D. Contribution 

There are two main problems investigated in the following 
sections, each leading to novel contributions; 

• We evaluate, for the hrst time, the statistical characteris¬ 
tics of the 7 d introduced in ([0. Due to its mathematical 
complexity, the exact pdf of 7 d is not derived but an 
approximation to the pdf is provided in this work. The 
approximated pdf is then used for calculating the outage 
probability and it is illustrated, by the simulations, that 
the approximated cdf leads to rather accurate results. 
Moreover, the exact cdf of 70 will be derived for the 
two asymptotic scenarios of full-correlation and no¬ 
correlation at the relay. Although this novel cdf is helpful 
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Fig. 1. Pdf of X with various number of A? paiameters. Monte Carlo 
simulations validate the correctness of the theoretical results. 


for outage analysis of the system, it is too complicated 
to be used for the analysis of the ergodic capacity. 

• In order to approximate the maximum achievable trans¬ 
mission rate, a very simple power allocation algorithm 
at the relay is introduced in this paper. As the optimal 
power allocation at moderate and high SNR is still an 
open problerrj^ for the purpose of comparison, exhaustive 
search over various discrete values of the rates is used as 
a benchmark. The values of the rates are obtained by 
allocating various amount of the power among different 
eigenmodes. According to the simulations, the proposed 
power allocation algorithm approximates the benchmark 
with insignificant difference. 


III. Statistical Analysis of Received SNR at 
Destination 


The SNR distribution in the relay is directly related to the 
ergodic capacity of the system evaluated in this paper. In 
order to maximize E{C(-)} in (j^, one should distribute the 
available relay power appropriately among different A^s in 


(14 1 so that E{C(-)} is maximized. Before we continue with 


the statistical characterization of 70 , two extreme scenarios 
are studied: we will investigate 70 when the relay does not 
experience any fading correlation and also when the antennas 
in the relay are fully correlated. These two scenarios, indeed, 
provide performance bounds, and so, the performance of a 
system with partial correlation will fall between the two 
bounds. 


^The optimal power allocation at low SNR was proposed in 


0 


A. Statistical Characteristics of^o Assuming Full Correlation 


Full correlation (FC) at the relay is equivalent to considering 
a system where all the elements of the S are unitj 0 i.e., 
S = Consequently, it is easy to see that Af = ttr and 

Xf = 0 for j > 2, and so, the random variable X in ([T^ can 
be written as 


Xfc 


AfAPATi 

1 + Af Af 


n\\^Xi 
1 + itrA^ Xi 


riR 


V 


1 + f^ 


(17) 


where, the random variable V = riRAf^Xi is exponentially dis¬ 
tributed with mean itrA^ . The following theorem introduces 
the cdf of 7 d for a fully correlated relay: 


Theorem 1. Assuming fully correlated antennas at the relay, 
the cdf of "fo is 


Pjo-rci^) = 1 - 2(nsu;)"^e-"*"' (18) 


ns—l 

xE 


771=0 


(XfusURw) / nsw ^ 

m\{ns — m—iy. y XfriR 


where Kj^(-) is the modified Bessel function of the second 
kind, i^-th order and w = 

"RTS 

Proof: See Appendix [A| for a detailed proof ■ 

The subscript “FC” in F,y^_p^{x) indicates the Full Corre¬ 
lation scenario. One can easily derive the pdf of 70 for the 
full correlation scenario by taking the derivative of ( [T^ with 
respect to x. 

In the next section, we study the statistical characteristics 
of the system for the non-correlated fading scenario. 


B. Statistical Characteristics of Jr, Assuming No Correlation 


No Correlation (NC) at the relay translates to S = 7. 
Such a scenario will occur when the relay antennas are placed 
sufficiently far apart and that the relay node is placed in a rich 
scattering environment. Assuming S = 7, it is straightforward 
to conclude that Af = Af = • • • = Ap = 1. On the 
other hand, since all Xj random variables follow the same 
distribution (exponential distribution with unit mean) and as 
all Xf values are equal to one, therefore all the Xf values 
should be assigned the same power, and so, let us assume 
Xf=X 2 =--- = X^^ = A®; consequently, the random vari¬ 


able X will be written as 


Xnc 


'^eq l^j=l 


V 

1 + V 


(19) 


where V follows the Erlang distribution with rate ^ and 

^cq 

shape nR. By substituting X^c in ( [1^ , the following expres¬ 
sion will be derived for the cdf of 70 for the non correlated 
scenario: 


Theorem 2. Assuming no correlation at the relay, the cdf of 


^We assume normalized correlation for simplicity of the notation. Gener¬ 
alization to a case with non-unity full correlation is straightforward. Similar 
normalization is assumed in Section [llI-B| tQO. 
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Fig. 2. Approximation of (^) as obtained in ^28) where = 4, 
A|^ = 1 and \G = \G = xG^ = 2. 



Fig. 3. Outage probability approximated by j28^ (ns oo) in comparison 
with the Monte Carlo simulations for various values of ns. 


Id is 


F'fo-.aix) = l-2{nswre-^^^ ( 20 ) 


riR — l ns — 1 

xEE 


m—O n—0 


nTi>^?qnsw) 

m\n\{ns — n — ly.w'^ m-n-i y 


with w = —. 

7s 

Proof: Following the same lines of the proof for ( [18) , 
one can easily prove ( |20) too. ■ 

Assuming single antennas at the source and the relay nodes, 
i.e., ns = riR = 1, the MIMO scenario of this paper reduces to 
a conventional single antenna relay system wherein the relay 
node has access to the variance of its channels; consequently, 
as expected, the two expressions in ( [TS) and ( |20) are identical 
according to 


V 7s7r V 7s7r 


) ( 21 ) 


Note that {x) for a single antenna scenario similar to ( |2T) 
has been reported in numerous papers including, e.g., | |22) . 

As discussed earlier, (and for arbitrary correlation 
should follow an expression that at the extreme case reduces 
to ( |20) and ( [TS) . 

On the other hand, statistical characteristics of X are also 
of essential importance for understanding the distribution 70- 
Assuming full correlation and no correlation at the relay, 
characterizing X was simple, however, characterizing X with 
arbitrary correlation is, actually, more complicated and will be 
derived in the next section. 


( fT^ is 


Fx{x) = <^ 1 - E n 

i—\ m = l 
m^j 


cc < 0 


Af+i < X < Xf 


Af < X 


where 


Ci = 


( 22 ) 


(23) 


A sketch of the proof for ( |22) is provided in Appendix To 
provide a better understanding of the statistical characteristics 
of the random variable X, the Fx{x) for k = 2 and 3 is 
provided in the following. Assuming k = 2, we have 


Fxix) = 


0 

1 - 
1 - 


-e ^ — 


X < 0 

0 < X < Af 
Af < X < Af 
X > Af 


(24) 


and for k = 3, Fx (x) is derived in (251 at the top of next 
page. One can easily derive fx{x) by taking the derivative of 
Fx{x) with respect to x, i.e., 

dFxix) 


fx{x) = 


dx 


(26) 


which, is a straightforward simple derivation practice. Fig. [T] 
illustrates the fx{x) assuming A^ = 1 for various values of 
AJ. The agreement between the theoretical and Monte Carlo 
simulations validates the correctness of the calculations. 


C. Statistical Characteristics of X for Arbitrary Correlation 

The statistical characteristics of the random variable X are 
not studied in the literature but will be derived in this paper. 
Considering that Af > Af > • • • > A^ > 0, the cdf of X in 


D. Statistical Characteristics of 7/j Assuming Infinite Anten¬ 
nas at the Source 

From the discussions provided in the previous section, it 
is clear that 7 d = 7sVX where Y follows an Erlang dis¬ 
tribution with rate and shape parameters equal to ns, i.e.. 
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ro 




a; < 0 

1 - 

C 2 C 3 pCiX 

C 1 C 3 

C 1 C 2 pC3X 

0 < a: < Af 

(C 2 -Ci)(c 3 -Ci) 

(C3-C2)(C1-C2) 

(C1-C3)(C2-C3) 

1- 

C 2 C 3 pCi x 

C 1 C 3 ^cox 


Af < a; < A 

(C2-Ci)(c3-Ci) 

(C3-C2)(C1-C2) 


1 - 

C 2 C 3 ^CtX 



Af < a: < A 

(C2-Ci)(c3-Ci) 



1 




a; > Af 


(25) 


fyiu) = and X follows a distribution as 

derived in ( |22| l and ( |2^ . In order to calculate the cdf of 7 d, 
one should calculate the following integral 

F^,{z) = P{7d < 4 = P{7si"2f < z} 

r°° z 

= / F{X <—} fY{y)dy. (27) 

Jo yis 


On the other hand, (271 does not lend itself easily to further 




Fx{-) 

Is 


(28) 


Proof: The random variable Y in the previous sections, 
e.g. in follows an Erlang-distribution; indeed, Y is the 
sum of ns exponential random variables, each with parameter 
ns (see ( |13a| l). Using the central limit theorem (see |24 
Ch. 7.4]), the random variable Y can be approximated by a 
Gaussian distribution with mean equal to 1 and variance 1/ns, 
i.e., approximately, Y ~ ^(1,1/ns). Assuming large ns (i.e., 
ns —oo), one can easily deduce U —> 1. By setting U = 1 in 
^n\, one can write F^^{z) = P{7d < z} = PjA" < —} and 


calculations, and to the best of our knowledge, the integral 
cannot be solved with the existing table of integrals (e.g., IH). 
However, clearly, F^^{z) is a multipartite function because 
X has a multipartite cdf; to the best of our knowledge, a 
multipartite F^^{z) in the context of AF cooperative systems 
has not been reported in the literature and so it is observed for 
the first time in this papei|^ The following theorem provides 
an approximate P/d(z) for the given system model; 

Theorem 3. For large ns, F^^{z) can be approximated by 


s^2^ is proved. 

Fig. and Fig. are intended to validate the precision of 
the approximation obtained in ( |28] l. In Fig. an illustration of 
Fj^ (z) is provided; clearly, Monte Carlo simulations approx¬ 
imate theoretical P/d (- 2 ) when ns is large. Moreover, con¬ 
sidering that the cdf in (|28|l corresponds to the outage 

probability for large ns (ideally tor ns 00 ), in Fig.we plot 
outage probability versus transmit power at the source node 
using ( [28] l and also using Monte Carlo simulations for various 
values of ns. It is clear that for large values of ns, the closed 
form expression for the outage probability approximates the 
Monte Carlo simulations with high accuracy. Nevertheless, for 

We use multipartite function to refer to a function that involves several 
distinct functions for different domains (e.g., see 124) and 125) ). Note that 
being multipartite is not considered to be advantage (or disadvantage) for 
the system but stressing on the novelty of the statistical characteristics of the 
SNR, in the context of AF cooperative systems, is meant to highlight the need 
for further investigation on the problem. 


smaller values of the ns, although the approximation is not 
accurate, it provides a reasonable approximation. 

As mentioned in |II-D[ one of the main objectives in this 
paper is to allocate available power in the relay according 
to (10 1 among different variables so that the ergodic 
capacity in (j^ is maximized. In fact, for plotting Fig. we 
arbitrarily assumed Af = A^ = 1; however, such a random 
power allocation to Af^ and A^ does not guarantee that the 
maximization problem in (j^ is solved. On the other hand, 
it was observed in ( |27l i that the Px(’) expression derived in 
( |22| is too complicated to lend itself to further mathematical 
calculations. Indeed, we do not know any closed form expres¬ 
sion for the objective function E{(7(-)} which can be used for 
calculating optimal A® values in (®. Furthermore, not only do 
we not know any analytical way for calculating optimal A^ 
values, we are not aware of any numerical method to calculate 
optimal Xf values. For the purpose of comparison, exhaustive 
search over various discrete values of the achievable rates is 
used as a benchmark. The rates for the exhaustive search 
are obtained by assigning various amount of the power to 
the eigenmodes according to the power constraint in 
moreover, the resolution of the exhaustive search is kept 
adequately small (0.1 dB) to ensure accurate approximation. 
Note that resolutions larger than 0.1 dB also provide accurate 
results, however, to make sure that no local maximum is 
missed, we use the the resolution of 0.1 dB throughout the 
paper when exhaustive search is provided for comparison. 

Although the proposed method in the next section is simple 
and straightforward, it will be revealed that the obtained values 
for A^s lead to the reasonable rates that are indistinguishable 
from the benchmark rates. Also, it will be revealed that the 
proposed method significantly reduces the computationally 
expensive calculations due to exhaustive search for finding 
optimal A^ values in real time practical communication sys¬ 
tems. 

IV. Two Antenna Relay 

For simplicity, as an initial step, let us assume a system 
with two antennas at the relay (i.e. tir = 2) where S = ^]; 

the parameter p indicates the correlation coefficient. As there 
are only two eigenvectors corresponding to S, the problem 
of optimal power allocation reduces to calculating the optimal 
values of A® and A^, given the power constraint in (10 1 . 

In fTT] Eq. 34], we derive a necessary and sufficient 
condition under which transmission only from the largest 
eigenvector (the eigenvector corresponding to Af) achieves 
capacity. For ease of reference, (m Eq. 34] is provided in the 
following lemma: 
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Fig. 4. Rate vs. transmit power at the relay. The rate values are obtained 
using the proposed power allocation algorithm, equal power allocation and the 
benchmark, ns = riR = 2, Nq = 1, Ps = 0 dB, and inter-antenna correlation 
p = 0.3. The exhaustive search is performed over discrete values of transmis¬ 
sion rates that were obtained using various values of and (step size 
0.1 dB) that fulfil 


Lemma: Transmission from the largest eigenvector achieves 
capacity if 


A 


with 


2 ^ 


(ai+PiA^)I?(Af,Pi)-aiE{ 


1 + Pl-^l 


} 


p gr 

-^1^1 1 + PiZi J 


(29) 


Note that ( [32l l and are obtained from the relay power 
constraint in ( [Tol l. ■ 

Remark: The motivation for assuming proportional power 
allocation in arises from the limit behaviour of the 

correlation coefficients. In the case where Af much larger than 
A 2 (i.e., Xf ^ A 2 ) clearly LER will be optimal transmission 
method and so Af^ ^ A^ must hold, and this is guaranteed 
by ( [33] l. On the other hand when Af and Af are only slightly 
different, both Af^ and A^ should be assigned relatively equal 
power; indeed, when Af = Af, t he fad ing is uncorrelated 
and so, as described in Subsection III-B Xf = Xf = A^, 


which again is guaranteed by (|33|).jNote that the conjecture 
will be validated in the following by simulations which show 
that the result is nearly identical to that with power allocation 
by exhaustive search. 

Fig. 0 illustrates the transmission rates of a cooperative 
system with ns = ur = 2 when the power allocation is carried 
out using the proposed algorithm. For comparison, maximum 
transmission rates corresponding to the benchmark are also 
illustrated. The figure clearly shows a good agreement between 
exhaustive search (i.e., the benchmark) and also the simple 
proposed algorithm. The difference between the proposed al¬ 
gorithm and the benchmark is, in fact, indistinguishable. Fig. 
shows the transmission rate assuming equal power allocation 
in the relay. Note that equal power transmission is equivalent 
to ignoring the knowledge of correlation at the relay. Clearly, 
the proposed algorithm significantly outperforms equal power 
transmission. 

With two antennas at the relay. Fig. shows that the 
proposed algorithm leads to excellent results. In the next 
section, the proposed algorithm is extended for a system with 
three antennas at the relay. 


22(Af,Pi) = 


= Xf{l + yX^Y)Xi 
1 1 


where Pi = 


Pr 


PiXf 
and ai = 


r(o, 


PiXf 


)e 


PlAj 


(30) 

(31) 


^ 2 Ps +-^0 


. Although 


111 


A^Ps+JVo ^ \]^Ps+No 

proves that FER is the optimal transmission method at low 

SNR, it does not discuss any method to distribute the available 
power at the relay among X^ variables when (29 1 does not 
hold. This problem is addressed in the rest of the paper. 


Proposition: Given Af, A 2 , Ps and Pr , if (29 1 holds, 
allocate the entire power in the relay only to Af^ as 


Af = 


No + PsXf 


(32) 


and set A^ = 0 (i.e., FER), otherwise, when (29 1 does not 
hold, we propose to allocate power per eigenvector propor¬ 
tionally to the strength of the eigenmodes, i.e., 

Af Af 


Af 


■ 


(33) 


Consequently, one can assume Af = Xfg 2 and Af = Afp 2 
with p 2 obtained from as 

Pr 


92 = 


Ps(AP + AP) + iVo(Af+ Af)' 


(34) 


V. Three Antenna Relay 

Assuming three antennas at the relay, it is clear that accord¬ 
ing to the values of Af, Af, Af, Pr and Ps, capacity optimal 
transmission can lead to three different scenarios: 

• Case 1: Transmission only via the largest eigenmode 
achieves capacity (i.e. transmission via the eigenvectors 
corresponding to Af ); or equivalently, FER is the 
capacity-optimal transmission method. This scenario will 
occur only when ( |29] l holds. In this case Af > 0 and 
Af = Af = 0. 

• Case 2: Transmission only via the two largest eigenmodes 
achieves capacity (i.e. transmission via the eigenvectors 
corresponding to Af and Af); or equivalently, 2-FER is 
the capacity optimal transmission method. In this case 
Af > 0, Af > 0 and Af = 0. 

• Case 3: Transmission via all three eigenmodes achieves 
the capacity; or equivalently, 3-FER is the capacity op¬ 
timal transmission method. In this case Af > 0, Af > 0 
and Af > 0. 

Note that ( |29l l specifies the FER-optimal region. In the 
following, we intend to specify necessary and sufficient con¬ 
ditions under which, assuming proportional power allocation, 
transmission only via the two largest eigenmodes in the relay 
(2-FER) approaches capacity. Note that proportional power 
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allocation is motivated by the precision of the algorithm 
introduced in Section IV We emphasis that the optimality 


of n-LER is conditioned on proportional power allocation and 
so it is sub-optimal, however, the results are acceptable when 
compared with the benchmarlj^ 

A. Conditional Optimality of 2-LER 

As discussed in Section [rv| when LER is not optimal, trans¬ 
mission via the two largest eigenmodes, with proportionally 
assigned power, approximates the benchmark with reasonable 
accuracy; therefore, by setting 

= >^^92 and = Af 52 (35) 

we aim to derive a necessary and sufficient condition under 
which, the maximization problem will be achieved by setting 
A® = 0 and assigning the available power in the relay, pro¬ 
portionally, to Af^ and A^. 


Considering that A^ > A^ > 


— •^nR’ 


one can 


easily conclude that if A^ > 0 leads to rate loss, then all 
the available power in the relay must be assigned only to 
Af and A^ and consequently A^ = 0 for j > 3. Let us assume 
that from the entire available power in the relay, e > 0 is 
assigned to A^ and, motivated by the results of the two 


antenna relay scenario in Section IV the rest of the power 
is proportionally distributed between Af^ and \9. It is easy 
to conclude from ^ -O) ^ q 

Cavils = 0) > Cavils = therefore, assigning e power to 
A^ will cause a rate loss. Now, let us assume that e —>■ 0, con¬ 
sequently, ^ < 0 is equivalent to ^ lAg->o — 0- 


dC 

Therefore, — ^ ^ |;^g_»q < 0 specifies a region in which 
assigning power to A® (and consequently A^ forj > 3) results 
in rate loss, and so one must transmit only via Af and A^ in 
this region. 

According to the power constraint in ( [TOl i, and assuming 
that A^ = Xfg 2 and A^ = X^g 2 , one can calculate the power 
assigned to Af^ and X 2 by calculating 52 in ([3^ as 


with 


P2 = 


Q !2 = 


g 2 — P 2 — 0 - 2 X 1 , 


Pr 


Ps(AP + AP) + iVo(Af+ Af) 
PsX? + No 


(36) 


(37) 


(38) 


Ps(AP + AP) + iVo(Af+ AP 
where the index of P 2 and a ;2 indicates that the 2-LER condi¬ 
tion is being considered. Substituting ( [T4| ) and ( [3^ in GD and 
considering that log(a/6) = log(a) — log(&), we obtain C{-) 
in at the top of the next page, where Wi = 1 -f 7 sVAf. 
The conditional optimality region of 2-LER corresponds to 
dC 

a region determined by —< 0’ that can be 
calculated by combining ( [39| ) and (|^ which is derived in ( |40| ) 
at the top of page with 


Z2 = APiEiATi -f ApiEzATz. 


S2t 


(42) 


^ Our conjecture is that the benchmark transmission rate obtained using 
exhaustive search is, virtually, equivalent to optimal transmission rate. 



Fig. 5. Blue area: Aj vs. A?. The area shows the LER optimal region. 

Red aiea: A? vs. AE The area illustrates the 2-LER optimal 
region for Aj = 0.5A^. Eor both the areas Pr = Pg = 10 dB. 


Note that the random variable X 3 at the second expectation 
operation on the right hand side of ( |40| l is independent of 
Z 2 and Y, hence, X 3 can be removed in the first expectation 
operation because £{ 263 } = 1. To the best of our knowledge, 
the first expectation and second expectations in ( |40l i cannot be 
further simplified; however, ( |40| can be further simplihed to 
( |4T] ) where 


V{Xl 2 ,P 2 ) = 


r( 0 :Ci) 

P 2 (Af-Ap 


r( 0 P 2 ) 

P2(Af - Af) 


(43) 


with Q = {P 2 Xf) Then, the conditional optimality region 
of 2-LER (i.e., ^ Ia^^-io — 0) can be obtained by some 

algebraic manipulation of ^4^ according to 


A^< 


(a 2 + PzAf )P(AP, P2) - a 2 nT+kz^} 


P 2 E{ 


l+XfjsY -i 

l + P 2-^2 ^ 


(44) 


Note that when ^ Ia^-io < 0 (i.e., not including 

equality), the necessary condition is also sufficient, and so, the 
strict “inequality” of ( |44l i specihes a necessary and sufficient 
condition, under which 2-LER is the optimal transmission 


method. However, when equality in ( |44l i holds (i.e., meaning 


that 


w 


IA 3 C- 


= 0), one should make sure that the 


optimum point is a maximum point; this can be done by show 


mg 

d^C, 


AA?) 


IA 3 °- 


< 0. In 111 


App. D], it is proved that 
ag_>.o < 0 is always valid. Eollowing the same lines 

< 0 


a(A?) 

C f A^) 

of proof, one can, similarly, prove that Ia^->o 

and so, the inequality in ( |44| i is a necessary and sufficient 
condition under which the 2-LER transmission is the optimaj^ 
transmission method. 


^Optimal in the sense that we assume proportional power allocation to 
the two largest eigenvectors. In the rest of the paper all n-LER transmissions 
are conditioned on proportional power allocation 
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C(-) — log + (P 2 — a 2 Xf)(^Xf^WiXi + A 2 ^T 4 ^ 2 ^ 2 ^^ ~ log + {P 2 — a 2 A^)(AfXi + AfX 2 ) + A^AfXa^. 


dCaviXf)^ 
dxG Ia?-o 


gr_a2/^| gf 

Xl + PoZoJ^ I 


. 1 + ^2-^2 

—E{-- 

P 2 M+P 2^2 


} + A^E{ 


Af(l + 7srAf)X3 

1 + P 2 Z 2 


£ , 1 + 7sFAj, ,02 


}-e{ 


}-( — 
1 +P 22’2 P 2 


Af X3 — a 2 {XfXi + Af'X2) 
1 + P2{x^ Xi + A^^2) 

Xf)V{Xl^,P 2 ) 


} 


(39) 


(40) 

(41) 


B. Near Optimal Power Allocation in a Relay with Three 
Antennas 


Similar to the algorithm of Section IV introduced for power 
allocation in a two antenna relay, a power allocation method 
will be introduced for the three cases {Case 1/2/3) discussed 
at the beginning of this section. 

Proposition: Allocate all available power in the relay 
to the largest eigenmode if ( [2^ holds (i.e., LER), otherwise 
check ( |44| l and allocate proportional power to the two largest 
eigenmodes if ( |44| ) holds (i.e., 2-LER); in the case when nei¬ 
ther (j2^ nor \AA) hold, then allocate the power proportionally 
to all three eigenmodes; i.e., Af^ = Afp 3 , A^ = A|'p 3 and 
Af = Afp 3 where 


Ps(AP + Ap + AP) + 7Vo(Af+ Af+ Af) 

The algorithm is summarized in Table ■ 

Fig. [^illustrates the transmission rates using the proposed 
algorithm. Comparison with the benchmark confirms that the 
proposed algorithm is effectively optimal. 


TABLE I 

Power Allocation in a Relay with tir = 3 


Step 1; Set Pi = 


~pr 


tVo+PsAf 


Step 2: Check the inequality in (29 1 
Step 3: If Step 2 is true 

Set Af = Pi, A^ = Ag^ = 0 and Quit. 
Step 4; Set P 2 from ( [J7| ) and 02 from ( |38| ) 

Step 5; Check the inequality in (44i 
Step 6 : If Step 5 is true 


Set Af = Xfg 2 , A® = Af 52 -, A^ = 0 and Quit. 


else 


Set Af = Ap 3 , X§ = Af 53 and A^ = Af 53 


VI. Proposed Power Allocation in a Relay with an 
Arbitrary Number of Antennas 

In this section, let us assume that the relay node is equipped 
with an arbitrary number of the antennas (say nf). As dis¬ 
cussed in earlier sections, in order to achieve capacity, the 
relay must assign an appropriate amount of its available power 
per eigenvector. Eollowing the same approach discussed in 
Section [V] one can assume k cases where, depending on the 
the system parameters (i.e., S, Ps and Pr), n-LER will be 
the capacity approaching transmission method; it means that 



Fig. 6. Rate vs. transmit power at the relay. The rate values are obtained 
using the proposed power allocation algorithm, equal power allocation and 
the benchmark, ng = 2, riR = 3, Ng = 1, Ps = 0 dB, and inter-antenna 
correlation pi 2 = 0.7, p 23 = 0.5 and pis = 0.2. The exhaustive seai'ch is 
performed over discrete values of transmission rates that wer e ob tained using 
various values of and (step size 0.1 dB) that fulfil jlO|. 


transmission via n (n < k) eigenvectors approaches capac¬ 
ity, and so, only n eigenvectors should be assigned power 
and the rest of the eigenvectors should be set to zero (i.e., 
X'f > 0 for j < n and A^ = 0 for j > n ). 

In Section IV we introduced the necessary and sufficient 
condition under which transmission via one eigenmode (largest 
eigenmode) achieves capacity; later on, in Section]^ a neces¬ 
sary and sufficient condition was derived, under which, trans¬ 
mission via the two largest eigenmodes approaches maximum 
transmission rate. One can extend the same concept and derive 
a necessary and sufficient condition under which transmission 
via n {n < k) eigenvectors will maximize the rate with the 
assumption of proportional power allocation. 


Analogous to Section [^ let us assume that from the 
available power in the relay, e > 0 is assigned to A^_|_i and 
the rest of the power is proportionally distributed between 
■ ■ ■ ^n- F is easy to conclude that 

is equivalent to the fact that assigning power to A^_,_i will 
result in rate loss, and consequently, since A® > A^ > • • • A^^^, 
the egienvectors corresponding to Xj for j > n + 1 should not 
be assigned power. 
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Theorem 4. The necessary and sufficient condition under 
which transmission from n largest eigenmodes (n-LER trans¬ 
mission) approaches capacity is: 


X 


n+l ^ 


where 


{an + PnXn+i)TX{Xf ... P„) Q!„E{ } 


P„E{- 




-} 


Pn = 


Pr 


and 


-A^O + Ps Xn+l 

" /\S ^n—2 

X'^'n 


(46) 

(47) 

(48) 


' An—2 

n(Af,... Pn) = Y. -r(0, (49) 

n (A^-A^) 


k = l 
k^m 


Proof: The proof is similar to that of Section V-A 


Note that ( |46l l is a general form of ( |29l ) and ( |44| . In fact, ( |46l ) 
will determine the eigenvectors that should be assigned power 
proportional to their correlation power. The power allocation 
algorithm is summarized in Table [H] The algorithm is indeed 
analogous to the water-filling algorithm: the purpose is to find 
the eigenmodes that should be assigned power for transmission 
and to discard the “weak” eigenmodes that result in rate-loss 
if assigned power. Note that the expressions derived in ( |29l l, 
( |44| l and ( |46| ) involve expectation operations that do not seem 
to lend themselves to calculation in closed-form; therefore, in 
practical implementation of the system one should implement 
them using numerical methods. Please see 13 wherein a 
numerical integration expression is derived. 

Fig. 0 illustrates the transmission rates for the proposed 
algorithm for a relay with four antennas and correlation 
coefficients as described in the caption. Clearly the proposed 
algorithm agrees with the rates obtained using exhaustive 
search. Moreover its superiority over equal power transmission 
is evident. The values chosen for inter-antenna correlation 
in the numerical simulations in Figs. and |7] are selected 
to be relatively high, and to reflect what might be expected 
in a relatively closely-spaced linear array. In these cases , 
the proposed algorithm demonstrates performance extremely 
close to the optimum. Note that when pij 0 (smaller 
inter-antenna correlation), the spatial correlation diminishes 
and so the equal power transmission will be the optimum 
method, which is indeed guaranteed by the proportional power 
allocation proposed in ( |3^ ; this is demonstrated by numerical 
simulation in Fig. Hence, since our algorithm is provably 
optimum at low correlation and is shown by simulation to be 
very close to optimum for a typical case of high correlation, 
it is at least a reasonable hypothesis that it is near optimum 
for all cases of practical interest. 


VII. Conclusion 

This paper studies the statistical characteristics of the re¬ 
ceived SNR at the destination in a MIMO relay network when 
the relay node experiences fading correlation. It is assumed 


TABLE II 

Power Allocation in a Relay with arbitrary ur 


Initiate n = 1 
while n < K 

step 1: Set P„ from ( |47| ) and a„ from ( |48l ) 
step 2: Check the inequality in ( |46l l 
step 3: If Step 2 is true 

Set for j <n 

Set A^ = 0 for j > n -f 1 and Quit while. 

else 

Set n ^ n -I- 1 and go to Step 1 
end (end of while when n > k) 



Fig. 7. Rate vs. transmit power at the relay. The rate values are obtained 
using the proposed power allocation algorithm, equal power allocation and 
the benchmark, ns = 2, nR = 4, No = 1, Ps = 0 dB, and inter-antenna 
correlation pi 2 = 0.7, pis = 0.5, pi4 = 0.3, p23 = 0.7, p24 = 0.5 and 
P 34 = 0.7. The exhaustive search is performed over discrete values of 
transmission rates that were obtained using various values of Ap (step size 
0.1 dB) that fulfil jT^. 


that the relay node has access only to the statistical CSI. In 
order to approach the ergodic capacity of the system, based 
on the available statistical channel knowledge at the relay, a 
new relay precoder design methodology is introduced. The 
proposed method is analogous to the water-filling algorithm; 
it searches for the largest eigenmodes that should be assigned 
power and discards the remaining eigenmodes. The simu¬ 
lations demonstrate good agreement between the proposed 
method and the benchmark which is obtained using exhaustive 
search. 
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Appendix A 

Proof of the SNR Distribution: Full Correlation 

By combining ( |T4l i and ( fTTl i, one can write 70 = 
and so, assuming w = we have 

= F{^<w)=F{V<y^) 

— w roo —-w 

I _ _ / (1 — 

Jo 

pOO -W 

1- / e^i’’">‘'“-”’/y(y)dy (50) 

Jo 

Note that V is exponentially distributed, and so, P(C < 
is non-zero only for F > 0; consequently, by substituting 
fyiv) in ( f5T| l and assuming the proper domain for the integral, 
we have 

00 

ns p ^ 

= 1 - J yns-i^-W^^r-^^dy (51) 


by changing variable according to f = y — w, assum- 

"S-l , 

ing + ^ “‘I ^PPly- 

m—0 


ing 123 3.471.9], F^^_p^{x) will be derived as 

F~i, 


7D-i.cH = l-2(nsw)”^e 


— nsw 


(52) 


{XfnsnRw) ('"+i)/2 


m—0 


m 


!(ns — 771—1)! 


nsw ^ 


and so, (|TF| is proved. 


Appendix B 

Statistics oe Random Variable X 
As a complete proof of ( |22l l is lengthy, we only prove the 
case of K = 2. Following the same approach, the extension 
of the proof to larger values of k (i.e., «; = 3,4, ••• ,n) is 
straightforward. Then by the rule of mathematical induction, 
it is easy to obtain the general expression in (|22li. 


Fx{x) for the Case of k = 2: 

The numerator and the denominator random variable X 
in include summation of random variables Xj which 
follow exponential distribution with unit mean. For simplicity. 


we assume k = 2 (i.e., X = 
of 




ArA 


X2 


l+A«AfXiA«A|=X2 


) and derive 


Fx{x) of ( |22l l. However, in order to simplify notation, in 
this section, let us substitute XfXj Xj consequently Xj 
is distributed exponentially with mean Xf. One can write 


Fx{x)=F{X <x) = 


\S2 V I \E2 V 
Ai A.I -h A 2 fV2 

1 + + A^-^2 


< x). 


(53) 


Applying basic algebraic manipulation, ( [5T| can be simplified 
according to 


p(a: <x) = p(Xi < 


+ Af (x — Af )2f2N 
Af(Af-a:) ' 


(54) 


Let us split the problem of deriving P(Ai < x) in over 
four different intervals: A) a; < 0, B) 0 < x < A2 C) Af < 
a: < Af and/)) Af < a; and derive P(X < a;) for each interval 
individually. 


A) X < 0.- 

Note that the random variables X, Xi and X 2 are non¬ 
negative random variables. Therefore, for a; < 0, we have 
¥{X <x)=0. 


B) 0 < a; < Af.- 

Assuming 0 < a; < A2 ; clearly, (Af — a:) at the denominator 
of ( |54| ) is positive. Therefore, since Xi > 0, the expression 
X -f A|‘(x — A 2 )Ar2 at the numerator of ( |5^ must be positive 
as well and, hence. 


2^2 < 


Ai^(x-Af) 

must hold. Therefore, (|54li can be simplified to 


(55) 


P(X < x) 


(56) 



X -f A2 (x — A|‘)x 2 
Af(Af -x) 


g 



fx2{^ 2 ) 


Considering that P(Ari < f) = 1 — e it is straightfor¬ 

ward to prove, from ( |56] l , that for 0 < x < Af, 

P(X < x) = 1--(57) 

C2 — Cl Cl — C2 

where Cj is defined in (|2^. 


C) A2 < X < Af .• 


It is easy to deduce from ( |5^ that for A2 < x < Af, the 
numerator and the denominator of ( |54l l is positive for every 
value of the random variable X^ and, hence. 


00 

p(a: <x) = j iP(^i < 


X -f A2 (x - Af )X2 ^ 


Af(Af-x) 


A? 


which leads to 

p(a: < x) = 1 - 


C2 


C2 - Cl 


dx2 (58) 


(59) 


D) X > Af.' 

For the given interval, clearly, the denominator of ( |54l l is 
negative while the numerator is positive; hence, one can easily 
conclude that assuming x > Af is equivalent to assuming 
Xi < 0. On the other hand, since Xi follows the exponential 
distribution, it is necessarily positive, and so, x > Af is an 
invalid assumption. Consequently, P(X < x) = 1 for x > Af. 
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